Measuring Qutrit-Qutrit Entanglement of Orbital Angular Momentum States 

of an Atomic Ensemble and a Photon 
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Three-dimensional entanglement of orbital angular momentum states of an atomic qutrit and a 
single photon qutrit has been observed. Their full state was reconstructed using quantum state 
tomography. The fidelity to the maximally entangled state of Schmidt rank 3 exceeds the threshold 
2/3. This result confirms that the density matrix cannot be decomposed into ensemble of pure 
states of Schmidt rank 1 or 2. That is, the Schmidt number of the density matrix must be equal to 
or greater than 3. 

PACS numbers: 03.65.Wj, 03.67.Mn, 32.80.-t, 42.50.Dv 
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An essential requirement toward practical quantum in- 
formation systems is the capability to generate entangled 
states between distant sites over quantum networks 
Inspired by a scheme to create long-lived entanglement in 
scalable quantum networks proposed by Duan et al. Q, 
various experimental results using atomic ensembles have 
been reported (such as 3]). Recently, d-dimensional 
quantum systems, or qudits, have been studied, and it 
has been pointed out that qudits are better adapted for 
certain purposes. As an example, they enable more effi- 
cient use of communication channels in quantum cryptog- 
raphy Q ■ Following the pioneering experiment on para- 
metric down-conversion ,5j, various protocols have been 
demonstrated using orbital angular momentum (OAM) 
states of photons (such as Q). Photons are a promis- 
ing carrier of quantum information. However, they are 
difficult to store for appreciable periods of time. The re- 
alization of massive qudits and the ability to characterize 
their entanglement are therefore critical for applications. 

Another recent landmark is the demonstration of the 
use of OAM to generate arbitrary superposition of atomic 
rotational states with the coherent transfer the OAM of 
light to atoms in Bose-Einstein condensate [7|, |8j. The 
experiment illustrates the potential of OAM as a tool to 
generate and to control the atomic qudits. 

Previous work Q demonstrated entanglement associ- 
ated with OAM in an ensemble of atoms and a photon. 
The atomic OAM state is linked to the spatial degree 
of freedom of collective atomic excitations, and, in the 
case of photons, the OAM state corresponds to Laguerre- 
Gaussian (LG) modes. This result suggests that atomic 
ensembles can be used as nodes of qudit-based quan- 
tum networks. However, previous observations were lim- 
ited to two-qubit entanglement. In this letter, we report 
higher-dimensionality of the entanglement of OAM states 



of an atomic ensemble and a photon, as confirmed by es- 
timating the Schmidt number [h| of the reconstructed 
two-qutrit (i.e., qudits with d = 3) density matrix. For 
pure states, the dimension of the range of the marginal 
state is called the Schmidt rank, which describes how 
many local levels are involved in the entanglement. Ex- 
tending this notion, the Schmidt number of a bipartite 
mixed state is defined to be the minimum Schmidt rank 
that must appear in any decomposition of the state as a 
mixture of pure states. For example, any decomposition 
of a mixed state with Schmidt number 3 must include a 
pure state of Schmidt rank 3 or greater. 

The LG modes constitute a complete basis for describ- 
ing the paraxial propagation of light 
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tensity and phase distributions of several LG modes and 
superpositions of them are shown in Fig.[T] An LG mode 
is characterized by its two indices p and m, and by the 
Gaussian beam waist wq. The integers p and m are the 
radial and azimuthal mode index, respectively, and the 
phase variation for a closed path around the optical axis 
is 2m7r. A single photon in the LG p ^ m mode carries a dis- 
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FIG. 1: The intensity (left panel of each pair) and phase (right 
panel) distribution of several LG modes and their superposi- 
tions. 
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crcte OAM of mh along its propagation direction. Here 
we define \L), \G), and \R) to be the single photon states 
with OAM of — ft, 0, and +ft, respectively. 

Consider an ensemble of atoms having a three-level 
structure \a), \b), and |c), as shown in Fig.^a). Initially 
all of the atoms are prepared in level \a). A classical write 
pulse tuned to the \a) — > |c) transition with proper detun- 
ing A is incident on the atomic ensemble. In this process, 
the |c) — > | b) transition is stimulated and a Stokes photon 
is generated. The write pulse is weak and its interaction 
time is short. Therefore the excitation probability of a 
Stokes photon into a specified mode is much less than 
unity per pulse. The collective atomic excitation retains 
the spatial distribution of the relative phase between the 
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FIG. 2: (Color online) (a) Energy levels of 87 Rb, and the 
associated laser frequencies, (b) Schematic of experimen- 
tal setup: SLM, spatial light modulator; SMF, single-mode 
fiber; SPCM, single photon counting module. Circularly 
polarized write and read pulses illuminate the S7 Rb MOT 
(Magnet-Optical Trap), and circularly polarized Stokes and 
anti-Stokes photons are selectively directed onto the SPCMs 
passing through quarter-wave plates and polarizers. The pair 
of SLM and SMF serves as a spatial mode filter, (c) Co- 
incidence rate and time-resolved coincidence (inset) between 
the Stokes and anti-Stokes photons of the LGo,o mode as a 
function of the time delay. 



write pulse and the Stokes photon. According to angular 
momentum conservation, the state of the Stokes photon 
and collective atomic excitation will be entangled. We 
use \l), \g), and |r) to denote the states of the collective 
atomic excitation with OAM of —ft, 0, and +ft, respec- 
tively. In the present work, our measurement is sensitive 
to only the three-dimensional photonic and atomic OAM 
states. 

When the write pulse carries zero OAM, The resultant 
atoms-photon state will be 

|0) = C L \L) \r) + C G \G) \g) + C R \R) \l) , 

where Cl, C'g, and Cr are the relative complex ampli- 
tudes. The amplitudes depend on the spatial shape of the 
region where the write pulse interacts with the atoms. 

The atoms-photon entanglement can be tested by map- 
ping the state of the atoms to that of an anti-Stokes 
photon by illuminating the atomic ensemble with a laser 
pulse (read pulse) resonant with the \b) — ► \c) transition. 
The efficiency of this transfer can be nearly unity be- 
cause it corresponds to the retrieval process of the atomic 
quantum memory based on electromagnetically induced 
transparency [131 ]. 

A schematic of the experimental setup is shown in 
Fig. Htb). An optically thick (optical depth of about 
5) cold atomic cloud is created using a magneto-optical 
trap (MOT) for 87 Rb. Levels |o) and \b) correspond to 
5Si/ 2 F = 1 and 2, respectively, and level |c) corresponds 
to 5Pi/ 2 F' = 2. One cycle in the experiment comprises 
a 6-ms loading period and a 4-ms measurement period. 
During the loading period, a gas of cold 87 Rb atoms is 
cooled and trapped, and they are optically pumped to the 
5S1/2-F = 1 level using a 50/is depumping pulse tuned to 



the 5S 1/2 F 



5P3/2F' 



2 transition. After the 



loading period, the magnetic field and the radiation re- 
sponsible for the cooling, trapping, and depumping are 
shut off. The vacuum cell is magnetically shielded using 
a single-layer permalloy. The coil jig is non-metallic, and 
thus eddy currents, which prolong the decay of the mag- 
netic field [§] , are suppressed. The residual magnetic field 
is about 1 mG during the measurement period. During 
that period, read and write pulses illuminate the atomic 
ensemble with a 400-ns repetition cycle. The 15-ns Gaus- 
sian write pulse is tuned to the \a) — > |c) transition with 
10-MHz detuning and comprises 4 x 10 4 photons. After a 
40-ns delay, a 200- ns rectangular read pulse of 300/iW in- 
tensity illuminates the MOT. As shown in Fig. [S^b) , the 
read and write pulses, whose spatial modes are cleaned up 
by passage through single- mode fibers SMF1 and SMF2, 
are counter-propagated along the same axis. The out- 
put beam from the fiber is focused into the MOT with 
a Gaussian beam waist of 400 /im. Similarly, the Stokes 
and anti-Stokes photons share a single spatial mode when 
their conversions are appropriately chosen using spatial 
light modulators SLM1 and SLM2 (Hamamatsu model 
X8267) . Their Gaussian beam waist at the center of the 
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MOT is adjusted to 275 /im. The angle between the axis 
of the Stokes or anti-Stokes photons and that of SMF1 
or SMF2 is ~ 3° in order to spatially separate the weak 
Stokes (anti-Stokes) photons from the strong write (read) 
pulses. The incident write and read pulses are circularly 
polarized, as are the Stokes and anti-Stokes photons. 

The Stokes and anti-Stokes photons coupled into 
SMF3 and SMF4 are directed onto the single-photon- 
counting modules SPCM1 and SPCM2 (Perkin-Elmer 
model SPCM-AQR-14). From the Stokes photon count- 
ing, the excitation probability is estimated to be 5 x 1CP 4 . 
Their outputs are then fed into the start and stop inputs 
of the time interval analyzer. The experimental results 
for the coincidence rate between the Stokes and anti- 
Stokes photons in an LGo,o mode are displayed versus 
time delay in Fig. [2c), from which the normalized cross- 

(2) 

intensity correlation is estimated to be g s ,&s = 74.6 ±7.4, 
confirming that the excitation probability of a Stokes 
photon into a specified mode is much less than unity for 
each pulse. The coincidence count rate was 5.2 s _1 . 

To determine the full state of the atoms and Stokes 
photon, two-qutrit state tomography [14j, was per- 
formed, where the density matrix was reconstructed from 
the set of 81 measurements represented by the operators 
jli <S> fij (with i, j — 0, 1, • • • ,8) and where p,k = \k)(k\. 
The ket | fc) for the Stokes photon was chosen from among 
{|L>, |G), \R), (\G) + \L))/V2, (\G) + \R))/V2, (\G) + 
i\L))/V2, (\G) -i\R))/V2, (\L) + |i?»/V2, (\L) + 
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FIG. 3: (Color online) Gaussian components of applied phase 
modulations T(x,y) to an incoming Gaussian beam of waist 
wo = 2.2 mm. (a) T(x,y) = e * ar s(— o+Kv-vo)) _ The 



left panel plots the simulation while the right panel shows 
the experimental result, (b) T(x,y) = e !arg(a: " 



T(x,y) = e *2 s s n ^-*o) In panels ( b ) and ( c ) 5 the dots are 
experimental results and the solid curves are obtained from 
the numerical simulation. 



i |i?))/\/2}, and for the collective atomic excitation from 
among |<7> , |r) , (\g)+\l))/V2, (\g) + \r))/V2, (\g)- 
i\l))/V2, (\g)+i\r))/V2A\l) + \r))/V2, {\l)-i\r))/^}. 
These measurements are implemented using SLMs and 
SMFs; the SLMs produce spatial phase modulation and 



the SMFs filter the LG , mode [16|. As reported in [17| 
in detail, arbitrary superpositions of an LGo.o and an 
LGq } ±\ mode can be converted into an LGo,o mode by 
application of the phase modulation T corresponding to 
the relative phase difference between the superposition 
mode and the LGo,o mode. As is clear from Fig. [TJ such 
conversions can be achieved by moving the singularity in 
the phase modulation to a particular location. Our re- 
flective SLMs have an active region of 768 px x 768 px. 
Even if the phase modulation is discrete, the fractional 
intensity diffracted into higher orders can be decreased by 
adding the blazed phase grating structure. The spatial 
period of the grating is 4 px ~ 100 /im with a diffrac- 
tion efficiency of 25 %. In order to check the SLMs, the 
Gaussian components of a beam diffracted with a spatial 
phase modulation of T{x, y) = e ^s(^x +t(y-y )) were 
measured. Here arg(z) is the argument of the complex 
number z. The results are plotted in Fig. [3ja) and (b), 
and are in good agreement with numerical calculations of 
the superposition modes. The location of the singularity 
that converts the superposition mode into a Gaussian can 
therefore be determined. At position (xo,j/o) = (0,0), 
where an incoming Gaussian mode is converted into an 
iGo,±i mode, the normalized intensity was measured to 
be 3 x 10~ 3 , indicating a high extinction ratio. Simi- 
larly, the superposition mode (LG ,-i + e' e LGo,+i)/\/2 
can be converted into an LGq^ mode by applying a dis- 
continuous phase modulation. The Gaussian components 
of the beam diffracted by an SLM with a spatial phase 



modulation of T(x, y) = e l 2 



if sgn(x— xo) 



was also measured, 



where sgn(x) is the sign of x. The experimental results 
are shown in Fig. [3jc), and are again in good agreement 
with numerical calculations. 

Fig. 2] shows the graphical representation of density 
matrix p exp reconstructed from the 81 coincidences. The 
typical coincidence rate was roughly 5 s _1 , and the data 
acquisition time of each measurement was 100 s. From 
the density matrix, the fidelity to a maximally entangled 
state F oxp = (MES|/5 oxp |MES) = 0.74 ± 0.02 was ob- 
tained. Here |MES) was chosen from the set of maximally 
entangled states (e M7r \L) \r) + |G) \g) + \R) 
so as to maximize the fidelity, where the values of a and 
(3 were 0.0197T and — 0.0587T, respectively. The error in 
i^cxp is calculated by using Monte-Carlo method from 
the statistical uncertainties in the coincidences. An op- 
timal witness operator of Schmidt number 3 in G 3 (£> G 3 
is given by W 3 = 1 - 3 |MES)(MES| /2 [H, resulting in 
Tv(W 3 p) < (MESIplMES) > 2/3. The experimen- 
tal result F cxp > 2/3 therefore confirms that the Schmidt 
number of the mixed state of the atomic ensemble and 
the photon is greater than or equal to 3. 
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FIG. 4: (Color online) Graphical representation of the den- 
sity matrix p cxp of a state as estimated by quantum state 
tomography from the experimentally obtained coincidences. 
The upper plot is the real part, and the lower plot is the 
imaginary part. 



the value 2000 expected from the measured excitation 
probability of 5 x 10~ 4 . Therefore, stray light appears 
to have strongly affected the photon statistics and 
decreased the fidelity. The imbalance of the diagonal 
elements may also affect the fidelity. However, even 
supposing that the Stokes photons are locally filtered 
so as to balance the relative amplitude, the fidelity is 
still expected to be 0.74. This result confirms that the 
imbalance is not the dominant factor decreasing the 
fidelity in our case. Note that the diagonal elements can 
be balanced by changing parameters such as the beam 
waist of the write pulse since the relative amplitude 
is dependent on the spatial shape of the effective 
interaction volume [19| . The other factor limiting the 
fidelity is the decoherence caused by the environmental 
noise, such as the Larmor precession of the ground-state 
Zeeman sublevels and the ballistic expansion of the 
atomic ensemble. 

In conclusion, higher-dimensionality of the entangle- 
ment of OAM states has been observed for an atomic 
ensemble and a photon by estimating the Schmidt num- 
ber of the reconstructed density matrix. The experiment 
described here enables one to communicate quantum in- 
formation encoded in the spatial deg rees of freedom of a 
photon and an atomic ensemble [201 ]. 
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The OAM measurements were achieved using mode 
conversion by SLMs and mode filtering by SMFs in this 
experiment. However, as frequently occurs in experi- 
ments utilizing a spatial phase modulation, the mea- 
surement bases cannot be realized completely accurately, 
resulting in unwanted radial and azimuthal components 
comprising up to 20%. While this systematic effect in- 
creases the error in F exp , the resultant fidelity is never- 
theless 0.74^°'^, which is larger than 2/3 even at the 
lowest error bound. 

The major diagonal elements of the recon- 
structed density matrix are (L\ (r|p cxp |L) |r) = 
0.25, (G\(g\p cxp \G)\g) = 0.37, and (R\(l\p exp \R)\l) = 
0.26. The summation of remaining elements 
1 - 0.25 - 0.37 - 0.26 = 0.12 suggests that compo- 
nents of non-zero total OAM, which may originate 
mainly from stray light, are one of the dominant factors 
limiting the fidelity. The normalized cross intensity 



correlation g : 



(2) 



74.6 ± 7.4 is significantly smaller than 
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